1  ^  AD-Aioe  043  analytic 
%  RESEARCH 
I  OCT  81 

unclassified  TASC-TR- 

SCIENCES  COUP  ftEAOINft  HA 
.‘2  If'iJJIi11  ,E?!I"ATI°N  *N0  CONTROL-OPTIMAL 

L 356-2 J  MVO  NOOOl" 

SAMPLE 

-79-C- 

F/G  1/3 

RA— ETC(U)  I 
0431  L 

NL  1 

1  11  1 

■ 

s 

L 

END 

l"'?2 

DTIC 

J 

MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  8URLAU  Of  STANDARDS  196.3  A, 


ONE  JACOB  WAY  READING,  MASSACHUSETTS  01867  (617)  944-6850 


81  12  02  036 


THE  ANALYTIC  SCIENCES  CORPORATION 


TR-1356-2 


RESEARCH  IN 

MULTIRATE  ESTIMATION  AND 
CONTROL-OPTIMAL  SAMPLE 
RATE  SELECTION 


8  October  1981 


Prepared  for: 

Office  of  Naval  Research 
Arlington,  VA. 


Prepared  by: 

Douglas  P.  Glasson 
Jeffery  Dowd 

Approved  by: 

James  D.  Shields 
E.  Wayne  Vinje 


THE  ANALYTIC  SCIENCES  CORPORATION 
One  Jacob  Way 

Reading,  Massachusetts  01867 


DlSTRIBUfiON  STA 


■■rrov'Xl  for  public 
iVlMrt ballon  Ontf.ii’tod 


UNCLASSIFIED 


security  classification  ok  this  aaoe  rv*n  o«» 


REPORT  DOCUMENTATION  PACE 


*C»0*T  HUMSCA 

TR-1356-2 


12.  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

).  RECIPIENT'S  CATALOG  NUMBER 


(PA3\ 


A.  TITLE  f«nW  SuMtllm) 

Research  in  Multirate  Estimation 
and  Control — Optimal  Sample  Rate 
Selection 


7.  AUTHO*r» 

Douglas  P.  Glasson 

i-  PERFORMING  ORGANIZATION  name  ano  aooress 

The  Analytic  Sciences  Corporation 
One  Jacob  Way 

Reading,  Massachusetts  01867 

"  CONTROLLING  OFFICE  NAMC  ANO  AOORESS 

Office  of  Naval  Research 
Arlington,  Virginia 

u  monitoring  agency  n  anC  *  ADORCisTif  SiimM  twmm  CrnOn 


».  TYPE  or  REPORT  A  PERiOO  COVCRCO 

5/1/79  -  5/1/80 

*■  PERFORMING  O AG.  REPORT  NUMBER 

TR-1356-2 

(.  CONTRACT  on  grant  numBEROJ 

N00014-79-C-0A31 


to.  AROGRAH  ELEMENT.  project,  TASK 
AREA  •  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

October  1981 

II.  NUMBER  OP  RACCS 


i  Controlling  Office;  1».  SECURITY  CLASS,  (mi  th»m  rmpmit) 

Unclassified 

l*e.  DECLASSIFICATION 'DOWNORADINO 
SCHEDULE 


I  »«  DISTRIBUTION  IT  ATBMCNT  (•<  thl*  Repeet) 


Approved  for  public  release;  distribution  unlimited. 


DISTRIBUTION  STATImINT  (ml  thm  etefrecf  on  from  In  Blmmt r  70.  it  dtiimrmml  irmm  *m*mrt) 


I  It.  SUPPLEMENTARY  NOTES 


Mi  KEY  WOK  OS  (Contlnum  mn  rmomrmo  e«#e  II  nmcmomnry  IBmntlty  Of  MecA  ntmmbmr) 


Digital  Control  Systems 
Flight  Control  Systems 


Multiple  Sample  Rate  Control 
Optimal  Control 


ACT  fCmnttnum  mn  revere#  aimm  it  nmmmmmnrr  mn4  tmmntltf  *r  MmJt  nmmimt) 


A  new  algorithmic  methodology  for  sample  rate  selection  in 
multirate  digital  control  design  is  developed.  The  methodology 
is  based  on  optimizing  the  random  disturbance  rejection 
properties  of  a  closed-loop  multirate  control  structure  under 
a  computational  constraint.  The  formulation  and  mathematical 
development  of  the  multirate  control  design  and  sample  rate 
selection  procedure  are  presented;  the  sample  rate  selection  _ 


DD  ,  'A??,  1473 


EDITION  ok  I  NOV  •(  It  OBSOLETE 


UNCLASSIFIED 


security  classification  ok  this  rase  i 


MCUMTV  CLASSIFICATION  W  THU  >*a«/TWf  Omm  truvl) 


l 


20.  ABSTRACT  (Continued) 


technique  and  properties  of  the  resulting  controllers  are  demon 
strated  through  aircraft  control  examples . 


UNCLASSIFIED 


THE  ANALYTIC  SCIENCES  CORPORATION 


PREFACE 


This  investigation  was  conducted  by  The 
Analytic  Sciences  Corporation,  Reading,  Massachusetts, 
from  1  May  1980  through  20  April  1981,  under  Contract 
N00014-79-0431  for  the  Office  of  Naval  Research, 
Arlington,  VA.  It  was  sponsored  by  the  Mathematics 
Program  headed  by  Dr.  Stuart  Brodsky.  Dr.  Charles 
Holland  was  the  cognizant  Navy  Scientific  Officer 
for  final  technical  review  and  public  release  of 
this  report. 

Mr.  Douglas  P.  Glasson  was  the  principal 
investigator  in  this  study.  He  was  supported  in 
this  study  by  the  technical  contributions  of  Mr. 
Jeffery  Dowd. 


iv 


V,  . 


4 


THE  ANALYTIC  SCIENCES  CORPORATION 


TABLE  OF  CONTENTS 


PREFACE 


Page 

No. 


iv 


1.  INTRODUCTION  1-1 

1.1  Background  1-1 

1.2  Summary  of  Results  1-2 

1.3  Organization  of  the  Report  1-4 

2.  OPTIMAL  SAMPLE  RATE  SELECTION  2-1 

2.1  Summary  of  the  Deterministic  Design  Technique  2-1 

2.2  Problem  Statement  2-7 

2.3  Mathematical  Formulation  of  the 

Optimization  Problem  2-10 

2.4  Optimization  Algorithm  2-13 

2.5  Prototype  Problem  2-17 

2.5.1  System  Description  2-17 

2.5.2  Optimization  Results  2-20 

2.5.3  Limitations  of  the  Gradient  Algorithm  2-23 

2.6  Chapter  Summary  2-27 

3.  EXAMPLE  SYSTEM  DESIGN  3-1 

3.1  Overview  3-1 

3.2  Mathematical  Formulation  3-1 

3.2.1  Aircraft/Disturbance  Dynamics  3-1 

3.2.2  Control  System  3-4 

3.2.3  Performance  Function  and 

Computational  Constraint  3-6 

3.3  Sample  Rate  Selection  Results  3-9 

3.4  Chapter  Summary  3-13 

4.  CONCLUSIONS  AND  RECOMMENDATIONS  4-1 

4.1  Conclusions  4-1 

4.2  Recommendations  4-2 

APPENDIX  A  COVARIANCE  AND  CONTROL  SOLUTION  SENSITIVITIES  A-l 
APPENDIX  B  MATRICES  FOR  THE  PROTOTYPE  PROBLEM  B-l 

APPENDIX  C  MATRICES  FOR  THE  EXAMPLE  SYSTEM  C-l 

REFERENCES  R-l 

DISTRIBUTION  LIST  DL-1 


| 


V. 

I 

t 

3 


1 

I 


5 

1 


v 


THE  ANALYTIC  SCIENCES  CORPORATION 


LIST  OF  FIGURES 


F igure 
No. 

Page 

No. 

2.1-1 

Multirate  Regulator  Structure 

2-3 

2.4-1 

Optimization  Algorithm 

2-14 

2.5-1 

Optimal  Sample  Rate  Selection  for  Space 

Shuttle  Example 

2-21 

2.5-2 

Performance  Function  in  C=0  Plane  for 

Prototype  Problem 

2-22 

2.5-3 

Payoff  Set  With  Reduced  Budget 

2-24 

2.5-4 

Payoff  Set  in  Nr=8  Plane 

2-24 

2.5-5 

Successive  Approximation  of  Performance 

2-26 

3.2-1 

Aircraft  Example  Models 

3-2 

3.2-2 

Flight  Condition  and  Performance  Function  for 
the  Example 

3-7 

3.3-1 

Performance  Function  for  Example 

3-10 

3.3-2 

Correlation  Coefficients  for  the  Example 

3-11 

3.3-3 

Performance  Function  at  Second  Flight  Condition 

3-12 

LIST  OF  TABLES 


Table 

Page 

No. 

No. 

2.2-1 

Problem  Statement 

2-8 

3.2-1 

Control  System  Variable  Weights 

3-6 

4 


j  v  .  ;  « 


i!*ri 


THE  ANALYTIC  SCIENCES  CORPORATION 


INTRODUCTION 


1 . 1  BACKGROUND 

The  practical  need  for  multiple  sample  rate  control 
systems  arises  from  the  finite  computational  capabilities  of 
digital  computers  used  in  control  implementation.  Although 
digital  computer  technology  continues  to  advance  significantly, 
particularly  in  the  areas  of  distributed  microprocessors  and 
very  large  scale  integration  (VLSI),  expanded  computational 
capabilities  are  inevitably  consumed  in  implementation.  In 
aircraft  applications,  for  example,  new  and  expanded  software 
requirements  for  such  functions  as  control,  navigation,  display, 
and  system  monitoring  tax  the  finite  storage  and  computational 
capabilities  of  the  flight  computers.  Accordingly,  the  flight 
control  system  designer  is  always  allocated  a  fixed  and  usually 
limiting  computational  capability  to  implement  a  control  design. 
In  addition,  accommodation  of  vehicle  flexibility  effects  (in 
the  form  of  instrument  output  filtering  or  active  structural 
control)  typically  demands  high  sample  rates  --  rates  that  may 
be  an  order  of  magnitude  higher  than  is  necessary  for  suitable 
control  of  the  rigid  body  modes  of  vehicle  motion.  Faced  with 
widely  varying  sample  rate  requirements  among  the  dynamic  modes 
of  the  vehicle,  a  multirate  control  structure  is  the  solution 
to  onboard  computational  limitations. 

In  the  first  phase  of  the  present  research  (reported 
in  Ref.  1),  a  systematic  methodology  for  designing  multirate 
control  systems  was  developed.  This  methodology,  which  is  based 
on  optimal  estimation  and  control  techniques,  obviates  dimension¬ 
ality-growth  problems  inherent  in  classically-based  techniques 
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(Refs.  2,  3,  and  4)  and  provides  a  step-by-step  procedure  for 
converting  a  desired  analog  control  design  to  an  equivalent 
multirate  design  without  approximation.  In  addition,  the  meth¬ 
odology  is  a  design,  as  opposed  to  analysis,  technique;  the 
iterative  trial-analysis-redesign  process  is  avoided  using  the 
new  technique. 

The  design  technique  developed  in  the  first  phase  of 
this  research  required  a  priori  specification  of  the  sample 
rates  for  the  system  controls.  In  the  present  report,  the  issue 
of  sample  rate  selection  to  optimize  system  performance  under 
a  computational  constraint  is  addressed.  A  general  mathematical 
formulation  of  this  optimization  problem  and  an  algorithm  struc¬ 
ture  to  solve  it  are  presented  and  illustrated  through  examples. 


1.2  SUMMARY  OF  RESULTS 


The  major  results  obtained  in  this  investigation  are: 


•  Basic  formulation  of  the  optimal  multirate 
estimator/regulator  design  problems 

•  Design-to-specification  procedure  for 
multirate  control  systems 

•  Extension  of  the  multirate  control  design 
formulation  to  proportional-plus-integral 
control  structures 

•  A  covariance  analysis  methodology  for 
periodic  control  structures 

•  Basic  formulation  of  the  optimal  sample 
rate  selection  problem 

•  A  general  algorithm  structure  to  solve 
the  optimal  sample  rate  selection  problem 

•  Investigation  of  optimal  sample  rate  solu¬ 
tions  and  closed-loop  system  properties 
through  examples 
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The  first  four  items  listed  were  developed  in  the  first  phase 
of  the  research  and  are  presented  in  detail  in  Ref.  1;  hence, 
these  items  are  only  summarized  in  the  present  report.  The 
remaining  three  accomplishments  are  the  primary  topics  of  the 
present  report. 

The  basic  formulation  of  the  optimal  sample  rate  selec¬ 
tion  problem  involves  the  development  of  system  performance 
metrics  and  computational  requirement  models  as  well  as  the 
combining  of  these  two  elements  into  a  Hamiltonian  function. 
Sufficient  conditions  for  a  stationary  point  of  the  Hamiltonian 
are  then  developed  with  due  regard  to  the  fact  that  integer¬ 
valued  independent  variables  (related  to  the  sample  rates  of 
slow  controls)  are  involved  in  the  formulation. 

A  gradient-like  minimum- seeking  algorithm  is  developed 
to  solve  the  optimization  problem.  This  algorithm  is  gradient¬ 
like  in  that  a  step-by-step  search  of  the  payoff  set  is  performed 
by  moving  in  steps  of  the  independent  variables  that  improve 
performance;  the  step  sizes,  however,  are  strictly  prescribed 
owing  to  the  fact  that  certain  of  the  independent  variables 
are  integer-valued. 

Two  example  applications  of  the  sample  rate  selection 
technique  are  presented.  The  first  example  is  the  lateral  dy¬ 
namics  of  the  space  shuttle  orbiter;  in  this  example,  the  basic 
properties  of  the  optimal  sample  rate  selection  methodology 
are  demonstrated  and  effects  of  restricted  computational  budget 
on  the  optimal  solution  and  resulting  closed-loop  system  proper¬ 
ties  are  evaluated.  The  second  example  is  optimization  of  the 
longitudinal  ride-qualities  of  the  F-14  aircraft  under  a  compu¬ 
tational  constraint.  In  the  F-lA  example,  realistic  models  of 
atmc  ’  *ri  turbulence,  which  constitutes  the  disturbance  to 
the  system,  are  used  and  physically-motivated  guidelines  for 
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selecting  the  appropriate  performance  metric  are  prescribed. 
Also,  variations  of  the  optimal  solution  with  flight  condition 
are  investigated  and  recommendations  are  made  for  practical 
control  implementation  based  on  the  observed  variations. 


1.3  ORGANIZATION  OF  THE  REPORT 

Sample  rate  selection  to  optimize  the  performance  of 
a  multiple  sample  rate  control  system  under  a  computational 
constraint  is  addressed  in  this  report.  An  overview  of  the 
deterministic  design  procedure  (developed  in  the  first  phase 
of  the  present  research)  is  presented  in  Chapter  2  in  addition 
to  the  problem  statement,  formulation,  and  solution  methodology 
of  the  optimal  sample  rate  selection  problem.  A  prototype  prob¬ 
lem,  based  on  a  lateral  control  design  for  the  space  shuttle 
is  also  presented  in  Chapter  2.  Application  of  the  new  sample 
rate  selection  technique  to  a  practical  aircraft  design  problem 
is  presented  in  Chapter  3.  Conclusions  and  recommendations 
for  future  research  are  presented  in  Chapter  4.  | 
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2. 


OPTIMAL  SAMPLE  RATE  SELECTION 


The  purpose  of  this  chapter  is  to  outline  the  mathemati¬ 
cal  formulation  and  solution  of  the  optimal  sample  rate  selection 
problem  and  to  illustrate  the  use  of  this  technique  in  a  prototype 
problem.  The  deterministic  part  of  control  system  design  (i.e., 
the  design- to-specif ication  procedure  to  obtain  desired  transient 
response  characteristics  from  the  multirate  control  system) 
was  addressed  in  detail  in  Ref.  1;  hence,  the  deterministic 
design  procedure  will  only  be  summarized  here.  The  problem 
statement  for  the  present  research  is  presented  in  Section  2.2; 
the  corresponding  mathematical  formulation,  derivation  of  opti¬ 
mality  conditions,  and  solution  methodology  is  covered  in  Sec¬ 
tion  2.3.  An  algorithm  to  solve  the  sample  rate  selection  prob¬ 
lem  is  described  in  Section  2.4.  A  prototype  problem,  based 
on  disturbance  rejection  optimization  of  the  space  shuttle  lateral 
controller,  is  solved  and  evaluated  in  Section  2.5.  The  chapter 
is  summarized  in  Section  2.6. 


2.1  SUMMARY  OF  THE  DETERMINISTIC  DESIGN  TECHNIQUE 

It  is  assumed  that  the  plant  dynamics  to  be  controlled 
are  described  by  the  linear-time-invariant  differential  equation: 

x ( t )  =  Fx(t)  +  Gu( t )  (2.1-1) 

where  x  is  the  state  vector  and  u  is  the  control  vector.  The 
control  u  is  to  be  sampled  and  held  at  time  instants  t^,  with 
all  elements  of  u  possibly  sampled  at  different  rates. 
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The  discrete-time  dynamics  of  the  continuous  plant  are 
described  by: 


-k+1  =  4>— k  +  rf“fk  +  rs-sk 


(2.1-2) 


where 


r  =  f  s  eFTG  dt 
Jo 


and  the  control  vector,  u,  has  been  partitioned  into  subsets 
U£  and  ug  to  signify  those  controls  computed  at  the  base  rate 
and  those  scheduled  at  lower  rates,  respectively. 


Development  of  the  multirate  regulator  structure  re¬ 
quires  augmentation  of  the  natural  plant  dynamics  by  holding 
circuit  states  (i.e.,  to  hold  the  slow  controls  between  updates); 
the  discrete- time  dynamics  of  this  augmented  system  (for  the 
case  of  two  sample  rates)  are  given  by: 


-k+1 

1 

4>  •  r 

1  s 
■ 

*k 

rf  !  rs6k,i£ 

| 

= 

1 

1 

1 

1 

-  -4 - 

1 

1  O 

_ 1 

-sk 

+ 

0  i  1  6k,i£ 

1 

(2.1-3) 


where  6k  is  the  Kronecker  delta  function. 


The  structure  of  a  multirate  optimal  regulator  for 
the  case  of  two  sample  rates  is  shown  in  Fig.  2.1-1.  As  the 
block  diagram  indicates,  one  control  channel,  uf,  is  updated 
at  a  fast  rate,  Tg  samples  per  second;  ug  is  computed  at  a 

slower  rate,  (£T  )"*  samples  per  second,  and  is  held  between 

s 


■i 
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Figure  2.1-1  Multirate  Regulator  Structure 
computations  by  a  holding  circuit.  Recomputation  of  u  is 

S 

accomplished  by  adding  an  increment,  v^,  to  the  holding  circuit 

on  cycles  when  k  =  i£ ;  on  all  other  cycles  is  set  to  zero. 

Another  result  of  the  present  design  procedure  is  that  the  slow 

control  is  crossfed  to  the  fast  control  channel.  The  purpose 

of  the  crossfeed  is  to  compensate  excitations  of  the  fast  modes 

of  the  plant  caused  by  u  on  cycles  between  slow  control  updates. 

s 

The  periodic  gains  --  Cf  ,  Cf  and  C  f  --  are  obtained 

rk  ri£  srk 

by  propagating  the  optimal  regulator  Riccatti  equation  from  in¬ 
finity  backwards  to  (periodic)  steady  state.  The  derived  peri¬ 
odic  regulator  gains  are: 


Ck  = 


Cfi£  :  Csfi£ 


ck  =  cf, 


k  =  i£ 


k  /  i£ 


(2.1-4) 


(2.1-5) 
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Note  that  the  dimensions  of  change  according  to  the  number 
of  controls  that  are  computed  on  cycle  k. 


The  closed-loop  dynamics  of  the  multirate  system  are 
determined  by  the  choice  of  the  quadratic  cost  functional  weight 
ing  matrices  used  in  the  optimal  regulator  problem;  i.e.,  with 
the  following  definition  of  the  cost  functional: 


t  <<4  si  )  Qk 

k=0  k  sk  k 

/x  k\ 

1  tv  l  m  m  A 

1  -"')  +  2(-k  -s,  )  Mk  -k 

sk 

^k  Rk  ^k} 

(2.1-6) 

The  choice  of  the  matrices,  Q^,  M^,  and  Rk>  ultimately  deter¬ 
mines  the  dynamics  of  the  closed-loop  system.  A  procedure  is 
described  in  Ref.  1  for  deriving  the  Q^,  M^,  and  R^,  matrices 
from  the  weighting  matrices  of  a  specified  continuous- time  regu¬ 
lator  design  (for  the  continuous  plant  dynamics  given  by  Eq. 
2.1-1).  This  procedure  is  summarized  here: 

(1)  A  continuous- time  regulator  is  designed 

to  meet  continuous-time  performance  speci¬ 
fications  . 

(•2)  Discrete-time  regulator  weighting  matrices 
for  the  single  rate  case  are  derived  from 
those  used  to  design  the  continuous- time 
regulator. 

(3)  The  discrete-time  periodic  weighting  ma¬ 
trices  used  to  solve  for  the  multirate 
gains  are  constructed  from  the  single 
rate  matrices.  The  multirate  gains  are 
derived  from  the  periodic  steady-state 
solution  of  the  discrete-time  Ricatti 
equation . 

The  first  step  involves  the  normal  iterative  design  process 
associated  with  designing  an  optimal  regulator;  i.e.,  trial 
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and  error  assignment  of  weighting  matrix  elements  to  yield  a 
regulator  that  meets  the  design  specifications.  The  relation¬ 
ships  presented  in  Ref.  10  are  then  used  to  compute  the  single¬ 
rate  (i.e.,  fast  rate)  discrete-time  weighting  matrices  from 
the  continuous- time  matrices.  Using  the  augmented  system  rep¬ 
resentation,  Eq.  2.1-3  and  the  relationships  in  Ref.  10,  expres¬ 
sions  are  derived  for  the  multirate  weighting  matrices;  inspec¬ 
tion  of  these  expressions  indicates  that  the  multirate  matrices 
can  be  constructed  from  the  single-rate  matrices  computed  in 
step  (2)  above.  An  implicit  benefit  of  this  three-step  approach 
is  that  the  deterministic  performance  of  the  continuous- time , 
single-rate  discrete- time ,  and  multirate  discrete-time  control¬ 
lers  can  be  compared  on  a  one-to-one  basis;  i.e.,  each  control 
design  is  based  on  optimizing  the  same  continuous- time  cost  func¬ 
tional  (Ref.  10). 


Extension  of  the  two- sample- rate  control  structure  de¬ 
scribed  here  to  the  multirate  case  (i.e.,  more  than  two  sample 
rates)  is  straight-forward.  Consider  organizing  the  plant  inputs 
into  two  groups:  those  computed  at  the  fastest  rate,  U£,  and 
those  scheduled  at  lower  rates,  ug.  Again  augmenting  the  plant 
dynamics  with  slow  control  holding  circuits  (which  now  are  oper¬ 
ating  at  multiple  rates)  the  discrete  time  dynamics  of  the  aug¬ 
mented  system  are  given  by: 


-k+1 
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where 


Ak  = 


6k,i1£1 


6 


k , ±2^2 
0 

•  0  0 


k,im,£m' 


where  m'  is  the  dimension  of  v  and 


,  a  15  Vi  =  k 

k*Vj  0;  i.£j  *k 


(2.1-7) 


(2.1-8) 


The  matrix  A k  is  a  generalization  of  the  Kronecker  delta  func¬ 
tions  used  in  the  two- sample- rate  case.  Each  diagonal  element 
of  nulls  the  input  to  a  particular  low  rate  control  holding 
circuit  on  cycles  between  new  computations  of  that  control. 

For  example,  if  us2  is  computed  every  two  cycles  (i.e.,  £ 2  -  2) 
then  the  second  diagonal  element  of  A^  will  be  zero  on  cycles 
when  k  is  not  a  multiple  of  two  (i.e.,  k  / 


Again,  the  deterministic  design  procedure  outlined  in 
this  section  maintains  consistency  of  transient  response  char¬ 
acteristics  in  proceeding  from  a  desirable  continuous- time  de¬ 
sign  to  a  multirate  digital  implementation.  In  the  following 
sections,  the  issue  of  optimizing  stochastic  error  rejection 
properties  under  computational  constraint  is  addressed.  Uni¬ 
formity  of  the  deterministic  properties  of  the  control  system, 
as  the  control  sampling  policy  is  varied,  is  implicit  in  the 
discussion  that  follows  as  a  result  of  using  the  design  proce- 
dure  described  here. 
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2 . 2  PROBLEM  STATEMENT 

Tne  issue  addressed  in  the  present  research  is  the  se¬ 
lection  of  sample  rates  for  a  multirate  digital  control  struc¬ 
ture  to  optimize  rejection  of  a  stochastic  disturbance  while 
accommodating  a  computational  constraint.  An  important  digital 
control  application  that  fits  the  general  problem  framework  de¬ 
scribed  here  is  optimization  of  digitally-controlled  aircraft 
ride  qualities  (performance)  while  flying  in  a  turbulent  atmos¬ 
phere  (disturbance)  under  computational  rate  limitations  of  the 
onboard  computer  (constraint).  In  subsequent  sections  of  this 
report,  the  formulation  and  solution  technique  described  here 
will  be  illustrated  in  typical  aircraft  applications. 

The  mathematical  problem  statement  is  presented  in 
Table  2.2-1.  The  plant  dynamics  are  described  by  a  linear-time- 
invariant  stochastic  differential  equation.  The  form  of  the 
plant  dynamics  shown  in  Table  2.2-1  includes  the  case  of  a  corre 
lated  process  (such  as  a  Markov  model  of  atmospheric  turbulence) 
disturbing  the  physical  plant  to  be  controlled;  in  this  case  the 
state  vector,  x,  includes  the  states  of  the  disturbance  process 
and  the  dynamics  matrix,  F,  includes  blocks  defining  the  distur¬ 
bance  dynamics  and  its  couplings  to  the  physical  plant. 

The  deterministic  control  design  is  based  on  optimizing 
an  infinite-time-horizon  quadratic  cost  functional,  J£ ,  (speci¬ 
fied  in  its  continuous -time  form  in  Table  2.2-1)  under  the  con¬ 
straint  of  piecewise  constant  controls  computed  at  multiple 
rates.  In  the  subsequent  formulation,  the  continuous-time  form 
of  the  quadratic  cost  functional  will  be  converted  to  its  equiv¬ 
alent  multirate  discrete- time  form  by  the  three-step  procedure 
described  in  Section  2.1;  hence,  uniformity  of  the  deterministic 
properties  of  the  control  design,  as  the  control  computation 
schedule  is  varied,  is  implicit  in  the  present  formulation. 
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TABLE  2.2-1 


PROBLEM  STATEMENT 


Given: 


Continuous-Time  Plant  Dynamics  (Linear-Time- Invariant 
Stochastic ) 

x  =  Fx  +  Gu  +  Aw  ;  w^  ~  N(0,o^) 

Control  Design  Cost  Functional  (Deterministic  Dynamics) 


J  =  i  f  {x^Qx  +  u^Ru}  dt 

c  z  Jq  -  - 


Error  Rejection  Functional: 


J  =  5  WkT 

Computational  Constraint: 
T  -j  m 


where 


f  =  fraction  of  computation  rate  capability  allocated 
c  to  control 

n  =  number  of  multiplications  per  control  channel 

t  =  execution  time  for  multiplication 
m 

Tg  =  base  sample  period 

N.  =  number  of  times  control  i  is  computed  per  control 
period 


Goal  of  the  Optimization: 

Determine  the  sample  period,  Tg ,  and  sample  policies,  N^, 

that  minimize  the  error  rejection  cost  functional  under 
the  computational  constraint 


V 


THE  ANALYTIC  SCIENCES  CORPORATION 


The  performance  to  be  optimized  is  expressed  by  the 
error  rejection  functional,  J.  In  the  form  shown  in  Table  2.2-1, 
the  performance  functional  represents  the  mean-squared  value  (in 
steady  state)  of  a  linear  combination  of  the  plant  and  distur¬ 
bance  states;  i.e.: 


=  Wk*k 


(2.2-1) 

(2.2-2) 


VkwI 


(2.2-3) 


Since  the  state  error  covariance  matrix,  P^,  of  the  closed-loop 
multirate  system  is  periodic  with  period  £ ,  a^  will  also  be 
periodic;  hence,  a  summation  over  £  cycles  and  normalization 
by  the  l/£  factor  are  included  in  Eq.  2.2-3.  The  elements  of 
the  weighting  array,  W^,  are  chosen  according  to  the  judgement 
of  the  designer  or  they  may  be  derived  on  a  physical  basis. 

In  an  aircraft  application,  for  example,  a^  may  be  a  physical 
quantity,  such  as  a  component  of  vehicle  acceleration;  the  ele¬ 
ments  of  can  then  be  derived  from  the  coefficients  of  the 
differential  equation  describing  that  acceleration. 

The  computational  constraint  is  a  simple  expression  of 
the  fact  that  a  finite  amount  of  computing  capability  is  avail¬ 
able  for  implementation  of  the  control.  The  first  term  of  the 
constraint  equation  is  the  computation  budget.  The  numerator 
of  the  first  term,  f cTg ,  is  the  portion  of  each  base  sample 
period  allocated  to  control  computations;  the  denominator,  ntm, 
is  an  approximation  of  the  computational  time  required  to  com¬ 
pute  one  control  channel  (i.e.,  element  of  u) ;  hence,  the  first 
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term  represents  a  specification  of  the  maximum  allowable  number 
of  control  computations  per  base  sample  period  (averaged  over 
£  cycles).  The  second  term  of  the  constraint  equation  is  the 
total  number  of  control  channel  computations  over  the  £  cycle 
control  schedule  (noim^iized  to  a  single  base  sample  period  by 
the  l/£  factor).  The  equality  expressed  by  the  constraint  equa¬ 
tion  stipulates  mat  the  total  computational  requirement  of  the 
control  implementation  should  equal  the  available  computational 
resources. 

The  gofrl  of  the  optimization  problem  is  to  determine 
the  values  of  the  independent  variables,  Tg  and  the  NL's,  that 
minimize  the  error  rejection  cost  functional,  J,  under  the  spe¬ 
cified  computational  constraint.  In  the  following  section,  the 
formulation  and  algorithmic  solution  of  this  optimization  prob¬ 
lem  are  described. 


2.3  MATHEMATICAL  FORMULATION  OF  THE  OPTIMIZATION  PROBLEM 


The  problem  statement  outlined  in  Section  2.2  can  be 
cast  into  a  formal  mathematical  framework  by  constructing  a 
Hamiltonian  function  that  couples  the  performance* measure  and 
the  computational  constraint.  The  form  of  the  Hamiltonian 
chosen  for  the  present  work  is  given  below: 

H  =  J  +  X  C  (2.3-1) 


i  * 

=  i  £ 

£  k=  1 


VkV 


+  X 


“c  nt 


m 


(2.3-2) 


The  Hamiltonian  is  constructed  by  adjoining  the  constraint  to 
the  performance  function  using  a  Lagrange  multiplier,  X.  Owing 
to  the  equality  constraint  that  C=0  at  the  stationary  point  of  H 
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H*  =  J*  (2.3-3) 

where  J  is  the  minimum  value  of  the  error  rejection  functional 
with  the  computational  constraint  satisfied. 

Finding  sufficient  conditions  for  a  minimum  of  the 
Hamiltonian  (and,  hence,  the  error  rejection  functional)  re¬ 
quires  consideration  of  the  fact  that  some  of  the  independent 
variables  of  the  problem  are  integer-valued  (i.e.,  the  NN's). 
The  optimization  strategy  in  this  case  will  be  to  use  the  usual 
zero  partial  derivative  condition  for  the  continuous  variables 
(T  and  X)  and  specify  that  the  Hamiltonian  should  be  nonde¬ 
creasing  with  respect  to  permissible  (integer)  changes  in  the 
' s  about  the  stationary  point.  Mathematically,  the  optimal¬ 
ity  conditions  are: 


1) 


3H 

3T 


i  * 

=  i  z 

1  i^i 


V, 


apk  „  T 
3TI  k 


+  A. 


(f«  "V ) = 


=  0  (2.3-4) 


2)  =  f  _JL 

BK  rc  nt. 


m 


m 


i  r 


ft  N‘ =  0 


(2.3-5) 


3)  H  Nondecreasing  with  respect  to  changes  in  N^'s 

At  this  point,  an  optimization  scheme  based  on  enumera 
tion  and  search  is  seen  to  be  a  possible  approach;  Eq.  2.3-5 
could  be  used  to  compute  all  permissible  combinations  of  the 
independent  variables,  then  all  combinations  could  be  tested 
to  determine  the  set  that  yields  the  best  error  rejection  (i.e. 
minimum  J).  In  the  present  work,  optimality  conditions  1  and  2 
are  combined  to  allow  more  flexibility  in  choosing  trial  sets 
of  the  independent  variables  and  the  search  over  the  trial  sets 
is  done  systematically.  The  development  of  this  approach  is 
outlined  here. 
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Using  optimality  condition  1,  given  by  Eq.  2.3-4,  a 
useful  interpretation  of  the  Lagrange  multiplier  and  Hamiltonian 
is  obtained.  Optimality  condition  1  can  be  written  as: 


9H  _  3J  .  3C  _  n 

3T  "  3T  A  3T  "  u 

s  s  s 


(2.3-6) 


.  _  3  J 

*  ~  ‘  3C 


(2.3-7) 


Using  Eq.  2.3-7,  the  Hamiltonian  can  be  rewritten: 


H  = J  -  tc  c 


(2.3-8) 


At  points  where  C  takes  on  nonzero  values; 


C  =  -  AC 


(2.3-9) 


where  AC  is  the  change  in  the  constraint,  affected  by  changing 
T_  alone,  required  to  satisfy  optimality  condition  2  (Eq.  2.3-5) 

s 

Hence : 


«  =  J  ♦  §£  *C 


(2.3-10) 


Equation  2.3-10  suggests  that  the  Hamiltonian  can  be  interpreted 
as  a  two-term  Taylor  series  approximation  of  the  error  rejection 
functional  with  the  N^'s  fixed  and  Tg  adjusted  to  satisfy  the 
computational  constraint.  This  Taylor  series  interpretation 
is  very  useful;  it  suggests  that  the  search  over  the  N^'s  can 
be  done  at  a  fixed  value  of  T  using  the  Hamiltonian  to  project 
the  performance  function  from  the  trial  point  (where  the  con¬ 
straint  may  not  be  satisfied)  to  the  appropriate  value  of  T  . 

An  optimization  algorithm  that  uses  this  interpretation  of  the 
Hamiltonian  is  described  in  the  following  section. 
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2.4  OPTIMIZATION  ALGORITHM 


An  optimization  algorithm  to  search  over  the  integer¬ 
valued  N^'s  is  shown  in  Fig.  2.4-1.  The  starting  point  is  a 
single-rate  controller  at  a  suitable  value  of  Tg  (i.e.,  all 
N.'s  equal  to  £  and  T  chosen  with  due  regard  to  the  sampling 

1  SO 

theorem).  The  performance  (Hamiltonian)  of  this  single  rate 
system  is  used  to  initialize  the  search  over  the  N^'s.  As 
shown  in  Fig.  2.4-1,  the  search  algorithm  has  an  outer-loop 
which  compares  the  values  of  the  Hamiltonian  obtained  from 
various  combinations  of  the  N^'s  (looking  for  a  minimum)  and 
an  inner- loop  that  generates  the  combinations  of  the  N^'s. 

The  inner-loop  operates  by  sorting  through  the  con¬ 
trols  one  at  a  time,  decrementing  N^  for  each  control  to  its 
next  lower  integer  value,  and  computing  the  Hamiltonian.  If 
the  Hamiltonian  for  the  decremented  is  lower  than  that  for 
the  previous  value  of  N^,  N^  is  set  to  the  decremented  value; 
if  not,  it  is  assumed  that  no  further  improvement  of  the  per¬ 
formance  can  be  obtained  by  changing  that  particular  (i.e., 
reducing  it)  and  its  value  is  locked  at  the  previous  value. 


The  outer-loop  iterates  the  operation  of  the  inner- 
loop  until  all  of  the  N^'s  are  either  locked  or  have  a  value 
of  1;  at  this  point,  no  further  improvement  of  the  performance 
is  possible  through  changing  any  of  the  N.'s  so  that  the  sta¬ 
tionary  point  has  been  obtained.  Finally,  the  N^'s  are  summed 
and  the  constraint  equation  (Eq.  2.3-5)  is  used  to  compute  the 
proper  value  of  Tg.  The  algorithm  shown  in  Fig.  2.4-1  is  a 
gradient-like  minimum- seeking  algorithm  with  prescribed  step 


sizes;  i.e.,  integer  changes  in  the  N^'s. 
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Computation  of  the  Lagrange  multiplier  used  in  the 
optimization  algorithm  requires  two  related  auxiliary 
computations : 

•  The  sensitivity  of  the  state  error  covari¬ 
ance  with  respect  to  Tg 

•  The  sensitivity  of  the  closed-loop  dynamics 
with  respect  to  T  . 

The  sensitivity  of  the  state  error  covariance  is  derived  by  dif¬ 
ferentiating  the  covariance  equation  with  respect  to  T  .  The 
covariance  equation,  which  has  a  periodic  solution  in  steady- 
state,  is  given  by: 


PW  =  *kVk  *  %  i  k=l,2,...,i 


(2.4-1) 


where 


6k  =  ♦  -  fkCk  =  <I-rkLk)6k 


(2.4-2) 


6k  =  *  •  rk  Rk  Mk 


(2.4-3) 


Differentiating  Eq.  2.4-1,  the  covariance  sensitivity  is  governed 
by: 


>Pk+1  _  3Pk  _T  96k  _t  -  aeJ  d% 

ir^karVar  PA  +  ekpk  ft  +  ft 

s  s  s  s  s 


with 


(I-rkLk>  ITT  *  W-  (I-rkLk>  6k 

s  s 
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Equation  2.4-4  is  a  periodic  Liapunov  equation  with  dynamics 
(0^)  corresponding  to  the  closed-loop  system  and  forcing  terms 


30, 


related  to  the  sensitivity  of  the  closed-loop  dynamics  (-=^) . 

Ts 

Once  a  (periodic)  solution  to  Eq.  2.4-4  has  been  obtained,  the 
Lagrange  multiplier  is  computed  from: 


A 


8P, 


W 


k  3T 


(2.4-6) 


The  sensitivity  of  the  closed-loop  dynamics,  required 

to  compute  the  covariance  sensitivity  (and,  hence,  the  Lagrange 

multiplier  and  Hamiltonian),  is  derived  by  differentiating  the 

control  solution  Ricatti  equation  with  respect  to  T  .  The 

s 

Ricatti  equation  is  given  by: 


K,.  = 


ekKk+l6k  +  Lk**kLk  +  Q: 


with 


L,.  = 


(rkVirk  *  V'1  rK+i 


(2.4-7) 


(2.4-8) 


Differentiating  Eq.  2.4-7  one  obtains  (after  algebra  covered 
in  Appendix  A) : 


— T  ^k+1  -  — t 

6k  ~3T  6k  +  6kKk+l 


h k  ark  \ 

\9Ts  '  9TS  Lk6kJ 


9rk  \T 

3Tg  Lkekj.  Kk+lGk 


+ 


3Q.  _  T  3R 

3T~  +  6kLk  3T"~  Lk6k 
s  s 


(2.4-9) 


2-16 


THE  ANALYTIC  SCIENCES  CORPORATION 


Equation  2.4-9  is  a  periodic  Liapunov  equation  with  dynamics 
corresponding  to  the  closed-loop  system  and  forcing  terms  re¬ 
lated  to  the  derivatives  of  various  plant  dynamics  matrices 
and  design  weighting  matrices.  The  (periodic)  solution  of  Eq. 
2.4-9  is  used  to  compute  the  sensitivity  of  the  closed-loop 


dynamics  matrix 
sensitivity . 


96. 

( — -) 
'8T  ' 


required  to  determine  the  error  covariance 


In  this  section  an  optimization  algorithm  for  select¬ 
ing  sample  rates  of  a  multirate  digital  control  system  has  been 
described.  In  the  following  section,  this  algorithm  is  applied 
to  a  prototype  problem. 


2 . 5  PROTOTYPE  PROBLEM 

In  this  section,  the  optimal  sample  rate  selection 
approach  developed  in  Sections  2.3  and  2.4  is  demonstrated  on 
a  prototype  problem:  a  lateral  controller  for  the  Space  Shuttle 
Orbiter.  A  description  of  the  prototype  is  given  in  Section 
2.5.1;  a  demonstration  of  the  optimization  technique  for  this 
example  is  presented  in  Section  2.5.2;  the  effects  of  restricted 
computational  budget  and  practical  design  considerations  that 
arise  from  the  observed  effects  are  discussed  in  Section  2.5.3. 

The  prototype  investigation  is  summarized  in  Section  2.5.4.  | 

2.5.1  System  Description  i 

I 

The  plant  to  be  controlled  is  the  lateral  dynamics  of 
the  Space  Shuttle  Orbiter  linearized  about  a  re-entry  flight 
condition : 

x  =  Fx  +  Gf6  +  G  6  +  Aw  (2.5-1) 

-  -  f  a  s  r  -  '  ' 
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XT  =  (v,  rb,  pb,  <t>)  (2.5-2) 

=  (sideslip  velocity,  yaw  rate,  roll  rate, 

roll  angle)  (2.5-3) 

=  aileron  deflection  (2.5-4) 

a 

6  =  rudder  deflection  (2.5-5) 

w  =  white  noise  turbulence  (2.5-6) 

E  { Aw( t+t )w( t )A }  =  Qc6(x)  (2.5-7) 


The  values  of  the  various  system  matrices:  F,  G^,  Gg , 
A  and  Qc  are  presented  in  Appendix  B.  It  is  worthwhile  to  note 
that  the  turbulence  model  used  in  this  example  forces  the  roll 
rate  dynamics  only;  i.e.,  Qc  is  of  the  form 


Q 


c 


0 

0 

0 


0 


0  0  0 

0  0  0 

0  qp  0 

0  0  0 


(2.5-8) 


The  deterministic  specifications  of  the  continuous¬ 
time  regulator  (from  which  the  multirate  designs  are  derived), 
in  terms  of  its  closed-loop  eigenvalues,  are: 


Roll-Spiral  Dynamics: 

s  =  -0.632  ±  i  0.913  (u»o  =  1.11  rad/sec,  £  =  0.569) 

(2.5-9) 


Dutch  Roll  Dynamics: 

s  =  -0.522  ±  i  0.159  (u j  =  0.546  rad/sec,  £  =  0.957) 

(2.5-10) 
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The  continuous-time  control  and  state  weighting  matrices,  R  and 
Q,  used  to  design  the  continuous-time  regulator  and  to  derive 
the  multirate  weighting  matrices  are  presented  in  Appendix  B. 

The  error  rejection  functional  used  in  the  present 
example  is  the  average  roll  angle  covariance  over  the  £  cycle 
control  period: 


i  s- 

=  \  E 

*  k=l 


VkV 


(2.5-11) 


with 

wk  =  (0,  0,  0,  1) 

J  =  j  E  p44  (2.5-12) 

k=l  k 

With  this  choice  of  error  rejection  functional,  the  sample  rate 
selection  technique  determines  the  sampling  policy  that  minimizes 
the  mean-square  deviation  from  wings- level  flight  caused  by  the 
turbulence  under  a  computational  constraint. 

The  control  period,  £,  and  computational  budget,  fc/ntm* 
for  this  example  are: 


£  =  8 


(2.5-13) 


£c_  . 
ntm 


=  6.25/sec 


(2.5-14) 


The  value  for  £  was  chosen  to  limit  potential  implementation 
complexity;  i.e.,  the  control  schedule  should  not  cover  a  large 
number  of  cycles.  The  computational  budget  parameters,  Eq. 
2.5-14,  were  chosen  to  illustrate  "interesting”  multirate  opti¬ 
mization  results.  It  should  be  mentioned  that  the  budget  chosen 
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is  much  tighter  than  that  used  in  the  operational  Shuttle  ve¬ 
hicle;  the  base  sample  rate  derived  in  the  subsequent  computa¬ 
tions  is  roughly  five  times  lower  than  that  of  the  operational 
Shuttle . 


In  this  example,  the  goal  of  the  optimization  is  to 
determine  the  base  sample  period,  T  ,  and  rudder  computation 
rate,  Nr>  to  minimize  J  (as  given  by  Eq.  2.5-12)  under  the 
computational  constraint  specified  by  Eq .  2.5-14.  The  aileron 
computation  rate,  N  ,  is  held  constant  at  8  (i.e.,  aileron  is 
computed  every  base  sample  period).  The  results  of  this  opti¬ 
mization  are  presented  in  the  following  subsection. 


2.5.2  Optimization  Results 

Demonstration  of  the  sample  rate  selection  technique 

described  in  Section  2.3  using  the  prototype  problem  defined 

in  Section  2.5.1  is  presented  here.  For  illustrative  purposes, 

a  survey  of  the  payoff  set  (i.e.,  values  of  the  error  rejection 

functional  over  a  range  of  values  of  the  independent  variables, 

T  and  N  )  is  presented  in  Fig.  2.5-1.  The  surface  of  the 
&  r 

payoff  set  is  smooth,  displaying  no  local  minimums  that  might 

"hang-up"  a  minimum  seeking  algorithm.  The  error  rejection 

function  is  seen  to  monotonically  increase  with  T  and  mono- 

s 

tonically  decrease  with  Nr  (recall  that  as  Nr  increases  the 
sample  frequency  for  the  rudder  increases);  the  behavior  of 
the  error  rejection  functional  makes  sense  in  light  of  similar 
behavior  in  single  rate  digital  systems  (Ref.  5). 

The  locus  of  the  constraint  function  (i.e.,  permissible 
combinations  of  Tg  and  Nr  under  the  computational  constaint) 

Although  it  is  shown  as  such  in  Fig.  2.5-1,  the  payoff  set  is 
not  a  continuous  surface.  The  payoff  set  is  comprised  of  a 
set  of  curves  for  fixed  integer  values  of  Nr- 
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Figure  2.5-1  Optimal  Sample  Rate  Selection 
for  Space  Shuttle  Example 

forms  a  plane  in  the  [Tg,  Nr>  J]  coordinate  system.  The  inter¬ 
section  of  the  constraint  plane  with  the  performance  "surface" 
is  the  payoff  set  of  interest;  the  optimum  error  rejection  under 
the  computational  constraint  is  located  on  this  intersection. 

The  method  of  stepping-along  the  payoff/constraint 
intersection  is  illustrated  in  Fig.  2.5-1.  As  described  in 
Section  2.3,  the  values  of  performance,  J,  along  the  payoff/ 
constraint  intersection  are  derived  (approximately)  by  comput¬ 
ing  the  Hamiltonian  along  the  Tg  =  0.18  sec  plane  and  interpret¬ 
ing  it  as  a  gradient-projection  of  the  performance  function 
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from  the  T  =  0.18  sec  plane  (with  N  fixed)  to.  the  C=0  plane. 
This  gradient  projection  of  the  performance  function  is  indi¬ 
cated  by  the  dotted  arrows  in  Fig.  2.5-1.  The  actual  minimum¬ 
seeking  is  achieved  by  step-by-step  decrementation  of  Nr 
(starting  with  a  single  rate  system,  Nr  =  8)  and  comparison  of 
the  Hamiltonian  function  for  each  successive  trial  value  of 
Nr>  In  proceeding  to  a  new  trial  value  of  Nr>  if  the  Hamil¬ 
tonian  is  larger  than  that  for  the  previous  value  of  Nr,  then 
the  optimum  has  been  overstepped,  and  the  stationary  point  is 
deduced  to  be  the  previous  value  of  N^. 

The  stationary  point  for  this  case  was  determined  to 
be  at  (Nf  =  2 ,  Tg  =  0.18)  as  indicated  in  Fig.  2.5-1.  The  opti¬ 
mum  and  its  neighboring  points  on  the  payoff/constraint  inter¬ 
section  are  shown  in  Fig.  2.5-2. 


R  72011 


NUMBER  OF  RUDDER  COMPUTATIONS.  Nr 

Figure  2.5-2  Performance  Function  in  C=0  Plane 
for  Prototype  Problem 

As  shown  in  Fig.  2.5-2,  there  is  a  gradual  improvement  of  per¬ 
formance  from  Nr  =  8  down  to  the  optimum  at  Nf  :  2  followed  by 
a  significant  increase  in  the  performance  function  (for  the 
small  decrease  in  Nf)  at  Nr  =  1.  Also  illustrated  in  Fig.  2.5-2 
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is  the  good  agreement  between  the  exact  values  of  the  perform¬ 
ance  function  and  those  derived  from  the  gradient  projection 
for  this  case.  In  the  following  section,  the  effects  of 
variations  in  computational  budget  (reduction  of  the  computa¬ 
tional  budget  in  particular)  on  the  optimal  solution  and  re¬ 
sulting  closed-loop  system  properties  is  investigated. 

2.5.3  Limitations  of  the  Gradient  Algorithm 

In  the  example  of  the  previous  section,  there  is  good 
agreement  of  the  exact  performance  and  the  approximate  perform¬ 
ance  via  the  two-term  Hamiltonian  expression.  When  the  budget 
is  tightened  (fc/ntffl  reduced),  the  Hamiltonian  may  not  provide 
a  good  projection  onto  the  constraint,  C=0 ,  plane. 

The  effects  of  a  restricted  computation  budget  are 
illustrated  in  Fig.  2.5-3.  The  computation  budget  has  been 
reduced  by  a  factor  of  four  over  the  previous  case  (Section 
2.5.2).  The  performance  function  in  the  C=0  plane  as  a  func¬ 
tion  of  rudder  computation  rate,  Nr,  is  shown  in  Fig.  2.5-3. 
There  is  a  sharp  optimum,  as  indicated  by  the  exact  perform¬ 
ance,  in  this  case.  The  performance  is  relatively  insensitive 
to  Nr  for  values  of  4  and  8;  with  less  than  4  (i.e.,  Nr  equal 
to  2  or  1)  there  is  a  rapid  degradation  of  performance.  The 
Lagrange  multiplier  approximation  is  inaccurate  for  Nr=4  and 
Nr=8,  causing  the  gradient-search  algorithm  to  incorrectly 
select  Nr=8  as  the  optimum  when  the  actual  optimum  is  at  Nr=4. 

Figure  2.5-4  shows  the  payoff  set  in  the  Nr=8  plane 
along  with  the  gradient  projection  onto  the  C=0  plane.  For 
the  example  illustrated  in  Fig.  2.5-4,  the  extrapolated  value 
of  J  at  C=0  is  a  poor  approximation  of  the  actual  performance 
due  to  the  curvature  of  the  payoff  set  between  the  computation 
point  (C=0.72)  and  the  constraint  plane  (C=0). 
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NUMBER  OF  RUOOER  COMPUTATIONS.  Nr 

Figure  2.5-3  Payoff  Set  With  Reduced  Budget 


VALUE  OF  CONSTRAINT  FUNCTION.  C 

Figure  2.5-4  Payoff  Set  in  Nr=8  Plane 
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There  are  many  possible  approaches  that  can  be  taken 
to  obtain  an  improved  projection  of  performance  onto  the  con¬ 
straint  plane.  One  such  approach  would  be  to  use  a  three- term 
"Taylor  Series"  approximation  of  the  form; 

JQ  =  J  +  AAC  +  j  AC2  (2-4-15) 

3  \ 

where  ^  is  derived  by  computing  A  at  two  values  of  Tg  and 
using  a  finite  difference  approximation. 

Another  approach  is  to  use  a  successive  approxi¬ 
mation  method  whereby  C  is  selected  according  to  a  modified 
Newton  procedure: 


Cn+1  =  ”Cn  + 


<VJn> 

PXr, 


(2.4-16) 


where  p  is  selected  to  reduce  the  extrapolation  interval  to 
the  C=0  plane.  For  the  case  of  p=2 ,  Cn+^  halves  the  interval 
to  the  C=0  plane;  thus,  determination  of  Hn+^  would  result  in 
an  improved  approximation  of  JQ. 

Use  of  the  successive  approximation  technique  with 
p=2  is  shown  in  Fig.  2.5-5.  An  improved  Hamiltonian  approxi¬ 
mation  is  obtained  with  the  interval  to  the  C=0  plane  reduced. 

The  successive  approximation  technique  requires  evaluation  of 
J  and  A  at  an  additional  point;  thus,  the  gradient  algorithm 
would  require  modification  to  include  this  additional  computation. 

At  this  point,  some  practical  design  issues  related 
to  optimal  sample  rate  selection  can  be  addressed: 

•  Reduced  budget  entails  a  larger  base 

sample  period  and,  hence,  overall  degrad¬ 
ation  of  closed- loop  performance 


a 


4r 
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VALUE  OF  CONSTRAINT  FUNCTION.  C 


Figure  2.5-5  Successive  Approximation  of  Performance 


Reduced  budget  may  result  in  a  controller 
having  high  sensitivity  to  parameter  varia¬ 
tions  (as  suggested  by  the  curvature  of 
the  payoff  set) 

The  designer  should  apply  judgement  in 
using  the  optimal  solution  when  high  sensi¬ 
tivity  (i.e.,  a  sharp  optimum)  is  indicated. 


In  comparing  Figs.  2.5-3  with  2.5-2,  an  overall  degradation  of 
performance  is  observed  for  the  restricted  budget  case  owing 
to  its  longer  base  sample  period.  The  sharp  optimum  evidenced 
in  Fig.  2.5-3  suggests  high  system  parameter  sensitivity;  the 
rapid  degradation  of  performance  for  reduced  rudder  sample  rates 
(reduced  N  )  suggests  that  at  the  optimum,  large  changes  in 
performance  may  result  from  relatively  small  changes  in  the 
design  model  parameters.  In  view  of  potential  sensitivity  of 


i  ' 
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the  optimal  performance,  the  designer  should  consider  imple¬ 
menting  a  suboptimal  sampling  policy  having  lower  sensitivity. 
For  example,  in  Fig.  2.5-3,  [Nr=8,  Tg=1.28  sec]  may  be  the  pre¬ 
ferred  sample  policy  for  implementation  since  it  should  be  less 
sensitive  to  system  parameter  changes. 


2 . 6  CHAPTER  SUMMARY 

In  this  chapter  the  problem  statement,  mathematical 
formulation,  and  solution  procedure  for  selecting  sample  rates 
for  multirate  control  systems  is  presented.  The  sample  rate 
selection  technique  optimizes  the  disturbance  rejection  proper¬ 
ties  of  the  multirate  controller  under  a  computational  con¬ 
straint;  an  equivalent  unconstrained  optimization  problem  is 
derived  by  a  Lagrange  multiplier  technique  and  a  gradient- like 
algorithm  to  solve  the  optimization  problem  is  described  and 
demonstrated  through  a  prototype  problem. 

The  gradient-like  optimization  algorithm  uses  a  Taylor 
series  interpretation  of  the  Hamiltonian  function  to  project 
the  error  rejection  performance  from  a  trial  point  (where  the 
computational  constraint  may  not  be  satisfied)  to  the  value  of 
the  base  sampling  period,  Tg,  that  satisfies  the  computational 
constraint.  The  algorithm  operates  by  systematically  stepping 
in  permissible  integer  values  of  the  slow  control  sample  rates 
to  find  the  minimum  of  the  performance  function  while  satisfy¬ 
ing  the  specified  computational  constraint. 

The  sample  rate  selection  technique  is  demonstrated 
through  a  prototype  problem.  With  a  reasonable  computation 
budget,  an  optimal  base  sample  period  and  slow  control  sample 
rate  are  derived;  error  rejection  performance  computed  by  the 
Lagrange  multiplier  approximation  and  the  exact  performance 
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are  in  close  agreement  for  this  case.  With  a  highly  restricted 
computational  budget,  a  strong  optimum  is  obtained  and  poor 
agreement  of  the  approximate  performance  and  exact  performance 
is  observed.  Modifications  to  the  optimization  algorithm  to 
improve  performance  projection  are  outlined  and  practical  con¬ 
siderations  in  using  the  derived  sampling  policy  in  implementa¬ 
tion  are  discussed. 

Again,  the  deterministic  design  procedure  (reported 
in  Ref.  1  and  summarized  in  Section  2.1)  along  with  the  sample 
rate  selection  technique  described  in  this  chapter  allow  the 
designer  to  optimize  the  disturbance-rejection  properties  of 
the  controller  under  a  computational  constraint  while  maintain¬ 
ing  desirable  transient  response  characteristics.  Specifically, 
the  flight  control  engineer  can  select  a  control  computation 
schedule  to  optimize  the  rejection  of  disturbances  (principally 
from  turbulence)  while  the  maintaining  specified  handling  quali¬ 
ties  and  satisfying  computational  constraints. 
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3.  EXAMPLE  SYSTEM  DESIGN 

3 . 1  OVERVIEW 

In  this  chapter  the  sample  rate  optimization  formula¬ 
tion  developed  in  Chapter  2  is  applied  to  a  practical  flight 
control  system  example.  A  dynamic  model  of  a  Navy  fighter  air¬ 
craft  with  a  multirate  control  system  structure  is  developed. 
With  an  appropriate  computational  constraint  and  performance 
function  defined,  the  optimal  sampling  policy  of  the  multirate 
system  is  selected.  The  performance  and  properties  of  the  opti¬ 
mum  system  are  investigated  as  a  function  of  flight  condition. 


3.2  MATHEMATICAL  FORMULATION 

3.2.1  Aircraft/Disturbance  Dynamics 

Figure  3.2-1  is  a  block  diagram  of  the  example  air¬ 
craft/disturbance/controller  model.  The  plant  to  be  controlled 
is  the  longitudinal  dynamics  of  the  F-14  aircraft  in  trimmed 
level  flight.  The  dynamics  of  the  vehicle  are  disturbed  by 
atmospheric  turbulence;  the  turbulence  model  used  in  the  present 
example  is  derived  from  the  Dryden  turbulence  spectra  (Ref.  6). 
The  aircraft  is  regulated  about  its  nominal  level-flight  condi¬ 
tion  by  a  multirate  proportional-plus-integral  controller. 

The  open-loop  continuous- time  aircraft/disturbance 
dynamics  are  given  by: 
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Figure  3.2-1  Aircraft  Example  Models 


x  =  Fx  +  Gu  +  Ad  w 


where , 


(3.2-1) 


T 

x  =  (6,  u,  q,  w,  wg,  x2) 

=  (pitch  attitude,  forward  velocity,  pitch  rate, 
vertical  velocity,  vertical  velocity  component 
of  turbulence,  auxiliary  turbulence  state) 

HT  =  <6e’  6mf> 

=  (elevator,  maneuver- flap) 

w  ~  N( 0 , 1 ) 


The  matrices  that  define  the  dynamic  model  of  Eq.  3.2-1  are: 


(3.2-2) 

(3.2-3) 
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F  =  vehicle  longitudinal  dynamics 


c  f  _  fi.  f  &  f 
±w  V  ±q  V  -q 


‘-p  p* 

.  0  _P_ 


v  ^  p  fq 


Ad  =  /3P  P 

72p  p 


(3.2-4) 


(3.2-5) 


(3.2-6) 


The  augmented  system  defined  by  Eqs.  3.2-1  to  3.2-6 
is  derived  by  coupling  the  turbulence  model  dynamics  to  the 
aircraft  dynamics.  The  dynamics  matrix,  F,  is  comprised  of 
blocks  related  to  the  aircraft  dynamics  (F),  the  disturbance 
dynamics  (Fd),  and  the  aerodynamic  coupling  between  the  dis¬ 
turbance  and  the  aircraft  (?c).  The  column  vectors,  and 
f  ,  are  derived  from  the  third  and  fourth  columns  of  F  and 
represent  the  aerodynamic  couplings  of  turbulence  velocities 
into  vehicle  motion.  The  scalar,  p,  is  related  to  the  tem¬ 
poral  bandwidth  of  the  turbulence  model  and  is  derived  from: 


b  -  5L 

P  L 

w 


(3.2-7) 


where 


V  =  aircraft  nominal  airspeed 
Lw  =  spatial  scale  of  Dryden  turbulence  model 


The  continuous-time  dynamics  given  by  Eq.  3.2-1  are 
converted  to  discrete-time  by  the  procedure  described  in  Sec¬ 
tion  2.1.  Numerical  values  of  the  various  system  matrices  are 
given  in  Appendix  C. 
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3.2.2  Control  System 

The  controller  used  in  this  example  is  a  multirate 
proportional-plus-integral  (PI)  controller  designed  to  the  same 
deterministic  specifications  as  the  F-14  control  system  studied 
in  Ref.  7.  The  gradient-search  algorithm  is  used  to  determine 
the  optimal  disturbance  rejection  of  the  closed-loop  system 
and  corresponding  sample  rates.  Since,  in  this  example,  the 
controller  is  to  maintain  trimmed  level  flight  as  accurately 
as  possible,  it  is  implemented  in  its  type  0  form;  the  feedback 
variable  transformation  to  obtain  a  pure  integration  in  the 
forward  control  path  (Ref.  1)  is  not  performed. 

The  controller  deterministic  specifications  are  em¬ 
bodied  in  the  elements  of  the  (continuous-time)  quadratic  cost- 
functional  matrices.  The  weighting  matrices  in  the  present 
example  are  generated  by  specifying  maximum  values  for  the 
various  state,  control,  control  rate,  and  output  variables  and 
using  the  following  formulae  to  relate  these  design  parameters 
to  the  states  and  controls  of  the  (augmented)  system: 


r  i 

i 

5 

Ax^  . 

max  l 

\ 

0 

V 

s 

1 

1 

gt 

Ax2 

max  iN 

L 0 

Au2 
max  i 

—  m 

V 

"\ 

1 

- 

Hx 

H 

u 

H„T 

Ay2 

•'max 

[f  5] 


(3.2-8) 
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(3.2-9) 


where 


and 


A  ■  M 

indicates  a  diagonal  matrix  with 


Axmax 

i 

“ 

maximum 

value 

of 

vehicle 

state  variable  i 

Aumax 

i 

= 

maximum 

value 

of 

control 

variable 

i 

Axmax 

i 

= 

maximum 

value 

of 

state  variable  i 

derivative 

Au 

max 

i 

= 

maximum 

value 

of 

control 

variable 

i  rate 

Ay 

Jmax 

i 

= 

maximum 

value 

of 

vehicle 

output  variable  i 

The  maximum  values  used  in  the  present  example  are 
listed  in  Table  3.2-1.  The  multirate  controller  weighting  ma¬ 
trices  are  computed  by  the  design-to-specification  procedure 
described  in  Section  2.1.  The  control  system  implementation 
is  given  by: 


(I-TsC2k  >Hk-i 


+  Ts  C3, 


(3.2-9) 


where  vg  k  is  the  slow  control  (maneuver  flap)  increment  computed 

at  an  integer  multiple  of  the  base  sample  period.  The  gains  -- 

C.  ,  C9  ,  Co  --  are  computed  from  the  solution  of  the  periodic 
1k  *k  Jk 

Ricatti  equation  corresponding  to  each  trial  control  computation 
schedule . 
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TABLE  3.2-1 

CONTROL  SYSTEM  VARIABLE  WEIGHTS 


MATRIX 

MATRIX  TYPE 

VARIABLE  NAME 

MAXIMUM  VALUE 

Q 

Vehicle  State 

Pitch  Angle,  6 

Axial  Velocity,  u 

Pitch  Rate, 

Vertical  Velocity,  w 

10  deg 

5  ft/sec 

20  deg/sec 

25  ft/sec 

Q 

Vehicle  Control 

Elevator  Deflection,  6 
’  e 

5  deg 

Maneuver  Flap  Deflection,  6  - 

mi 

00 

Q 

State  Derivative 

Pitch  Angle,  0 

Axial  Acceleration,  U 

Pitch  Acceleration,  qb 

Vertical  Acceleration,  w 

0D 

25  ft/sec^ 

00 

00 

Q 

Output 

Normal  Acceleration,  a 

Maneuver  Flap  Deflection,  6  , 

mi 

3  ft/sec^ 

5  deg 

R 

Control  Rates 

Elevator  Rate 

Maneuver  Flap  Rate 

3  deg/sec 

6  deg/sec 

3 

2.3  Performance 

Function  and  Computational 

Constraint 

The  steady-state  variance  of  the  vehicle  normal  accel¬ 
eration  at  the  vehicle  center  of  gravity,  an>  in  response  to 
vertical  turbulence  is  chosen  as  the  performance  function;  the 
physical  environment  corresponding  to  this  choice  of  perform¬ 
ance  index  is  depicted  in  Fig.  3.2-2.  Normal  acceleration  is 
an  appropriate  measure  of  vehicle/control  system  performance 
on  physical  grounds  in  that  it  is  closely  related  to  ride  qual¬ 
ity  standards  (see  Refs.  8  and  9,  for  example).  Maintaining 
vertical  accelerations  at  acceptable  low-levels  is  an  important 
control  function  with  respect  to  preventing  pilot  fatigue  and, 
possibly,  motion  sickness. 
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Figure  3.2-2  Flight  Condition  and  Performance 
Function  for  the  Example 


Mathematically  the  error  rejection  performance  function 
is  defined  by: 


J  =  E{a/} 


(3.2-12) 


with 


where 


VkV 


(3.2-13) 


a 


n 


=  Vq  -  w 


(3.2-14) 


P.  =  steady-state  covariance  matrix  of  the 
closed-loop  system 


W.  =  weighting  array 


V, 
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£  =  control  schedule  period 


V  =  nominal  airspeed 


and  E [ * ]  is  the  expectation  operator. 


The  elements  of  the  weighting  array,  W^,  are  derived 
from  Eq.  3.2-14  using  the  vehicle  state  equation  for  w,  i.e.,: 


w  =  lx  + 


X 

i  =  W,  --- 

n  k  u 

L  ~  J 


=  V(0, 0,1, 0,0, 0,0,0)  -  (?  g) 


(3.2-15) 


(3.2-16) 


(3.2-17) 


where  1  and  g  are  the  row  elements  of  the  F  and  G  matrices 
(Eq.  3.2-1)  which  constitute  the  differential  equation  for  w. 

For  this  example,  the  computational  constraint  (Table 
2.2-1)  reduces  to: 


fcTs  1 

-  F  (N  +  N  c)  =  0 
nt  £  e  mf 

m 


(3.2-18) 


where 


number  of  elevator  control  computations  per 
£  cycle  period 


number  of  maneuver- flap  computations  per 
£  cycle  period 


The  control  period  and  computation  budget  are: 

£  =  8 


(3.2-19) 
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Y~  =  20.0/sec 
m 


(3.2-20) 


The  control  period  value  was  chosen  in  consideration  of  impler 
mentation  complexity;  i.e.,  a  control  period  longer  than  8  cycles 
would  entail  involved  logic  in  the  software  to  "keep  track  of" 
the  control  schedule.  The  computation  budget  was  chosen  to 
achieve  a  50  percent  reduction  of  computation  over  a  single¬ 
rate  controller  operating  at  a  20  Hz  sample  rate. 


3.3  SAMPLE  RATE  SELECTION  RESULTS 

The  results  of  the  performance  optimization  for  the 
F-1A  example  are  presented  in  this  section.  Several  flight 
conditions  are  chosen  to  determine  system  performance  and  sample 
rate  sensitivities  to  variations  in  flight  condition.  Although 
an  exhaustive  survey  of  the  entire  F-1A  flight  envelope  was 
not  performed,  it  is  possible  to  derive  some  useful  conclusions 
from  the  flight  conditions  selected. 

The  error  rejection  performance  of  the  aircraft  in 
trimmed  level  flight  at  an  altitude  of  20,000  ft  and  an  air¬ 
speed  of  999  ft/sec  is  shown  in  Fig.  3.3-1.  A  strong  optimum 
is  evidenced  compared  to  other  flight  conditions  surveyed;  the 
optimal  sampling  policy  is  Nmj  =  A  and  Tg  =  0.075  sec.  There 
is  a  slight  reduction  in  mean  square  normal  acceleration  as 
Nraf  is  decreased  from  8  to  A;  further  reduction  in  the  number 
of  maneuver- flap  computations  results  in  significant  degrada¬ 
tion  of  performance;  thus,  the  maneuver- flap  gust  alleviation 
potential  rapidly  degrades  at  low  sample  rates  (Nm£<A). 
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R  72184 


NUMBER  OF  MANEUVER  -  FLAP  COMPUTATIONS.  Nmf 

Figure  3.3-1  Performance  Function  for  Example 

The  correlation  coefficients  for  maneuver- flap  and 
elevator  against  vertical  velocity,  and  Pw/ge>  and  against 

the  vertical  velocity  turbulence  component,  P5f/Wg  and  Pfte/Wg> 
for  different  maneuver- flap  computation  rates  are  shown  in  Fig.  i§ 

3.3-2(a),  (b).  The  cross  correlations  rapidly  drop  off  past 
the  optimal  (Nm£<4);  this  indicates  that  at  low  sample  rates 
the  flight  control  system  loses  its  ability  to  track  distur¬ 
bance  inputs  due  to  the  large  time  period  between  maneuver- 
flap  updates. 

The  error  rejection  performance  in  trimmed  level  flight 
at  an  altitude  of  20,000  ft  and  an  airspeed  of  400  ft/sec  is 
shown  in  Fig.  3.3-3.  For  this  flight  condition,  the  performance 
is  insensitive  to  sampling  policy  (i.e.,  the  system  “xhibits  a 
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very  weak  optimum.  The  optimum  occurs  at  N  ,  :  2  and  T  =  0.0625 
sec .  The  performance  insensitivity  to  sampling  policy  illustrated 
in  Fig.  3.3-3  indicates  that  a  potential  reduction  in  total  com¬ 
putation  time,  due  to  reduced  maneuver-flap  sampling  rate,  can 
result  from  a  multirate  design;  hence,  multirate  design  offers 
a  distinct  computational  advantage  over  single  rate  for  this  case. 
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OPTIMAL  SAMPLING  POLICY 
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Figure  3.3-3  Performance  Function  at 
Second  Flight  Condition 

A  practical  design  procedure  to  derive  a  sampling  policy 
to  cover  a  complete  flight  coindition  envelope  can  be  conjectured 
at  this  point: 


•  Perform  an  exhaustive  survey  of  optimal 

sampling  policies  as  a  function  of  flight 
condition  over  the  vehicle  operational 
envelope 
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•  Determine  the  subset  of  flight  conditions 
that  exhibit  strong  optimums 

•  Choose  a  single  control  schedule  that 
has  minimal  off-optimal  performance  de¬ 
gradation  over  the  entire  flight  envelope. 


In  the  present  example,  choice  of  the  optimum  schedule 
for  20,000  ft  at  999  ft/sec  would  not  result  in  significant 
performance  degradation  for  the  other  flight  condition  surveyed 
(e.g.,  for  h  =  20,000  ft  at  v  =  400  ft/sec,  a  suboptimal  design 
would  be  acceptable  due  to  low  sensitivity) 


3 . 4  CHAPTER  SUMMARY 

In  this  chapter,  the  sample  rate  selection  procedure 
has  been  applied  to  an  example  aircraft:  An  F-14  that  is  "flown" 
through  a  vertical  atmospheric  disturbance.  Steady-state  variance 
of  the  closed-loop  aircraft  normal  acceleration  at  the  center-of- 
gravity  is  chosen  as  the  performance  index.  The  sample  rate 
optimization  problem  for  the  case  of  a  fixed  computation  budget 
is  solved  to  determine  the  optimum  choice  of  sample  period, 

Tg,  and  maneuver- flap  computation  rate,  Nm£. 

Variations  of  the  optimal  solution  as  a  function  of 
flight  condition  were  examined.  It  was  found  that  about  strong 
optimums  there  is  a  rapid  drop-off  of  correlation  between  control 
action  and  disturbance  input  as  the  slow  control  sample  rate 
is  reduced.  Other  flight  conditions  surveyed  show  performance 
that  is  insensitive  to  sampling  policy.  A  practical  design 
procedure  would  be  to  perform  an  exhaustive  survey  of  the  entire 
flight  envelope,  to  look  for  flight  conditions  having  strong 
optimums,  and  to  choose  a  control  schedule  based  on  the  strong 
optimums . 
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4 .  CONCLUSIONS  AND  RECOMMENDATIONS 


The  net  result  of  the  present  research  (including  that 
reported  in  Ref.  1)  is  a  general  technique  for  digital  flight 
control  design.  Using  the  design- to-specif ication  procedure, 
the  designer  is  able  to  synthesize  multirate  digital  controllers 
having  uniform  handling  qualities  (i.e.,  deterministic  response) 
which  are  independent  of  the  control  schedule  chosen  (within 
the  limitations  of  the  sampling  theorem).  Using  the  sample 
rate  selection  procedure  outlined  in  the  present  report,  the 
disturbance  rejection  properties  of  the  controller  are  opti¬ 
mized  under  a  quantitative  computational  constraint  (i.e.,  com¬ 
puter  specifications  and  budget  are  included  as  design  param¬ 
eters).  As  a  result  of  these  two  (automated)  procedures,  the 
designer  can  conveniently  structure  and  design  a  control  system 
entirely  in  continuous- time ;  conversion  to  a  multirate  digital 
implementation  and  selection  of  control  sample  rates  to  accommo¬ 
date  computational  constraints  is  then  achieved  automatically. 


4.1  CONCLUSIONS 


The  research  reported  in  this  document  provides  the 
following  significant  results: 


•  Formulation  of  the  Optimal  Sample  Rate 
Selection  Problem  -  The  optimal  sample 
rate  selection  problem  is  formulated  as 
a  constrained  optimization  problem  and 
converted  to  an  equivalent  unconstrained 
problem  by  Lagrange  multiplier  techniques. 
General  formulae  for  the  error  rejection 
performance  functional  and  the  computa-' 
tional  constraint  are  presented;  sufficient 


4-1 


V 

•  ■■  v.  f  ;>V,,X,U'JW«,  '  1 vm  . 


THE  ANALYTIC  SCIENCES  CORPORATION 


conditions  for  a  stationary  point  of  the 
corresponding  Hamiltonian  function  are 
derived. 

•  Optimization  Algorithm  -  A  gradient-like 
algorithm  to  solve  the  optimal  sample 
rate  selection  problem  is  derived.  This 
algorithm  performs  a  step-by-step  optimiza¬ 
tion  of  the  cost  functional  moving  in 
prescribed,  integer-valued  steps  of  the 
control  sample  rates.  Limitations  of 

the  algorithm  for  highly  constrained  compu¬ 
tation  budgets  are  investigated;  alternative 
measures  to  overcome  these  limitations 
are  discussed. 

•  Example  Application  -  The  sample  rate 
selection  technique  is  applied  to  a  prac¬ 
tical  design  example:  a  flight  control 
system  for  the  F-14  aircraft.  A  physically- 
motivated  performance  measure,  rms  normal 
acceleration,  is  selected  and  quantitatively 
related  to  the  general  error  rejection 
functional.  Optimal  control  schedules 

are  derived  for  a  limited  set  of  flight 
conditions;  practical  considerations  for 
control  schedule  selection  are  made  based 
on  these  results. 


4 . 2  RECOMMENDATIONS 


It  is  recommended  that  the  following  studies  related 
to  the  present  investigation  be  undertaken. 


•  Investigate  the  effects  of  control  schedule 
on  control  system  sensitivity  to  uncertain 
plant  parameters.  Determine  appropriate 
parameter  sensitivity  metrics.  Formulate 
and  derive  a  solution  methodology  to  opti¬ 
mize  the  sensitivity  metric  under  a  computa¬ 
tional  constraint.  Investigate  the  properties 
of  the  optimal  solution  through  a  design 
example . 
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•  Apply  the  techniques  developed  in  the 
present  research  to  a  practical  flight 
control  system.  Aeroelastic  vehicle  con¬ 
trol  design  would  be  a  particularly  attrac¬ 
tive  test  application. 

•  Investigate  the  application  of  multirate 
control/sample  rate  selection  techniques 
to  other  areas  of  control;  e.g.  ,  process 
control  scheduling,  multi- target  tracking, 
remote-to-central  packet  switching  schedule, 
and  attention  allocation/dynamic  decision 
modeling  in  human  operator  applications. 
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APPENDIX  A 

COVARIANCE  AND  CONTROL  SOLUTION  SENSITIVITIES 

Derivations  of  the  covariance  and  control  sensitivity 
equations  (Eqs.  2.4-4  and  2.4-9,  respectively)  are  presented 
here.  The  covariance  sensitivity  equation  is  derived  in 
Section  A.l;  the  control  solution  sensitivity  is  presented  in 
Section  A. 2. 


A.l  COVARIANCE  SENSITIVITY 

As  discussed  in  Section  2.4,  direct  dif ferentation  of 
the  covariance  equation,  Eq.  2.4-1,  results  in 


3P 

3T 


k+1  . 

s 


—  ®^k  T  _*p  _  ^®k 

ek  FT  ek  +  FT  Pkek  +  ekPk  FT  *  FT 


(2.4-4) 


with 


<1“rkLk> 


(1‘rkLk> 


} 


(2.4-5) 


It  is  necessary  to  evaluate  the  first  factor  of  the  second 
term  in  Eq.  2.4-5: 


(I-rkLk> 


8T 


3  'I-rk<rkKk+lr^k>'lrkTKk*ll 


§f;  <Itrk“kV,W'1 


(A. 1-1) 

(A. 1-2) 
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Equation  A. 1-2  is  a  more  convenient  form  for  this  derivation 
and  was  obtained  using  the  matrix  identity: 


[I-YT(ZYT+X)_1Z]  =  (1+YTX-1Z)‘1 


(A. 1-3) 


with 


x  =  Rk 

y  =  z  =  rkK*{* 


Carrying  out  the  dif ferentation  in  Eq.  A. 1-2  one  obtains: 

sr  “♦ViV’W'1  *  -(I+rkRk,rIKk.1)"  5T  l(Itrk“k‘riVi>l 

s  s 


(A. 1-4) 


-  -(,-rkLk>  ar  la'ViVVi'W-w' 

(A. 1-5) 


The  matrix  identity  of  Eq.  A. 1-3  was  used  again  (in  reverse)  to 
derive  Eq.  A. 1-5  from  A. 1-4.  Carrying  out  the  differentiation 
of  the  inner  factor  of  Eq.  A. 1-5  one  obtains: 


w  (I-W  =  -(I’rkL 
8 


f|"^k  '-1  T  *  ^k  ~-l  '-1  il 

kJ  Rk  Pk+rkRk  8^  ’  PkIRk  ^  Rk 


A-1  ®^k+ll  T 

*  rkRk  rk  iF  |  «-w 
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36.  3Qd 

The  remaining  derivatives,  and  gy, 

s  s 

straightforward  differentiation: 


are  derived  by 


(A. 1-7) 


(A. 1-8) 


l 

£ 


I 


v 


(A. 1-9) 

(A. 1-10) 
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Differentiating  Eq.  A. 2-2,  one  obtains 

3Kk  T  3Kk+l  30k  T  30k  3^k 

3T  ‘  Bk  3T  0k  3T  Kk+l0k  °k  Rk+1  3T  3T 

s  s  s  s  s 

(A. 2-4) 

The  partial  derivative  of  is  derived  by  differen¬ 

tiating  Eq.  A. 2-3. 


3Kk*i 


(I-rkLk> 


(I-rkLk> 


Kk*i  sr  rk<rJIWVV 

s 


3Fk  T  ark 

Kk+1  3T  Lk  '  Kk+lLk  3T 
s  s 


‘lrTK  I 

r kKk+ll 


Lkrk  ^Sr1  +  Kk+lrk  FT  <rkKk+lrk+^k)  lrkKk+l 

s  s 


(A. 2-6) 


■&l  <rX.irk*,t1'1  * 


<rkKk+irktRk> 


l8Ts 


T  3rk 

Kk+irk+rkKk+l  3T“ 


.  rT  aKk+l  r  8Rk 

1  k  3T  1  k  3T 
s  s 


}  <rkK 


k+irk*Rk> 


(A. 2-7) 


so  that 
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*ei  -  (I.r  L ,  .  K  ^  t  .  lt  *s  K  .  L  rt  3V. 


3T  3T 
s  s 


k  k  k+1  3T  k  k  3T  "k+1  "k  k  3T 
s  s  s 


'I3rk 

C  1.3Ts  Kk+irk 


T  a^k  T  a^k+l 
+  rkKk+l  3T~  +  Pk  _3T  rk  + 

6  S 


»kl 

aTs|Lk 


(A. 2-8) 


or 


aKk+l  9Kk+1  T  T  ark 

3T  "  (I‘rkLk)  3T  (I"rkLk)  '  (I“Lkrk)Kk+l  3T  Lk 
s  s  s 


lk|af=  Vi  <-’-r 


kLk>  *  ir]1 


(A. 2-9) 


Substituting  Eq.  A. 2-9  into  Eq .  A. 2-4  one  obtains 
the  final  result: 
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APPENDIX  B 

MATRICES  FOR  THE  PROTOTYPE  PROBLEM 


The  matrices  which  describe  the  vehicle  dynamics,  dis¬ 
turbance  dynamics,  and  control  design  of  the  example  of  Sec- 
tion  2.5  are  presented  here  : 


Continuous  System  Dynamics: 


F 


-0 . 3970E-01 
-0.2970E  00 
-0.2590E  01 
0.0 


-0.1000E  01 
-0 . 4120E-01 
0.1140E  00 
0.0 


0.2094E 
0 . 5000E-02 
-0.2940E  00 
0.1000E  01 


0 . 2207E-01 
0.0 
0.0 
0.0 


(B-l ) 


Fast  Control  Effectiveness  (Aileron): 


G 


f 


-0 . 4981E-03 
0.1184E-01 
0.2390E  01 
0.4369E-04 


(B-2) 


Slow  Control  Effectiveness  (Rudder): 


G 


s 


0.4125E-02 
-0.1649E  00 
0.8405E  00 
-0 . 4828E-04 


(B-3) 


*The  numerical  values  presented  here  are  given  in  standard  "E 
format";  the  notation  EsNN  indicates  a  factor  of  10s  .  For 

example,  0.1000E-02  is  equivalent  to  0.001. 
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Process  Noise  Covariance: 


Qd  = 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.2500E  00 

0.0 

0.0 

0.0 

0.0 

0.0 

(B-4) 


Continuous-Time  State  Weighting  Matrix: 

0.1000E-01  0.0  0.0 

0.0  0 . 2500E-02  0.0 

0.0  0.0  0.2500E-02 

0.0  0.0  0.0 


Q  = 


0.0 
0.0 
0.0 

0 .  A000E-01 

(B-5) 


Continuous -Time  Control  Weighting  Matrix: 

n 


R  = 


0.2500E  00 

0.0 


0.0 

0.2500E  00 


(B-6) 


$ 


a 
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APPENDIX  C 

MATRICES  FOR  THE  EXAMPLE  SYSTEM 


The  matrices  which  describe  the  vehicle  dynamics,  dis¬ 
turbance  dynamics,  and  control  design  of  the  example  of  Sec¬ 
tion  3  are  presented  here  (Flight  Condition:  V  =  999  ft/sec; 
h  =  20,000  ft). 


I 

C.l  AIRCRAFT/DISTURBANCE  DESCRIPTION 
(Eqs.  3.2-2  THROUGH  3.2-5) 

Continuous-Time  Aircraft/Control  Dynamics: 


f  = 


0.0 

-0.5615E  00 

0.0 

-0. 2050E-02 


0.0 

-0 . 1069E-01 
-0 . 1681E-02 
-0 . 5946E-01 


0.1000E  01 
-0 . 8830E-01 
-0.1141E  01 
0.1754E  02 


0.0 

0.4301E-01 
-0.6618E  00 
-0.1375E  01 


(C.l-l) 


Disturbance  Dynamics: 


-0.5709E  00 

0.0 


0.5709E  00 
-0.5709E  00 


(C.  1-2) 


*The  numerical  values  presented  here  are  given  in  standard  "E 
format";  the  notation  EsNN  indicates  a  factor  of  10s  .  For 

example,  0.1000E-02  is  equivalent  to  0.001. 
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Coupling  Dynamics: 


f 


c 


0.0 

0 . 4590E-01 
-0.6244E  00 
-0.1949E  01 


0.0 

-0 . 2891E-02 
-0.3736E  01 
0.5744E  00 


Continuous-Time  Control  Effectiveness: 


G  = 


0.0 

0.2358E  00 
-0.1805E  02 
-0.2554E  01 


0.0 

-0.3895E  00 
-0.1135E  02 
-0.4691E  01 


Process  Noise  Influence: 

0.0 

0.3050E  01 
-0 . 3739E-01 

A  .  = 

a  -0 . 3302E-03 
-0 . 4266E-02 
0 . 6560E-01 


(C. 1-3) 


(C. 1-4) 


(C. 1-5) 


C. 2  CONTROL  DESIGN 

l 


Continuous-time  Augmented  State  (vehicle  state  and 
control)  Weighting  Matrix: 


C-2 
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0. 1050E-01 

0.2315E-04 

0.5575E-04 

0. 2745E-03 

0.3699E-03 

0. 1418E-02 

0.2315E-04 

0.4039E-01 

-0.6823E-03 

0.9083E-02 

0.1687E-01 

0.3100E-01 

0.5575E-04 

-0.6823E-03 

0.3703E-02 

-0. 1583E-01 

-0.2940E-01 

-0.5389E-01 

0.2745E-03 

0.9083E-02 

-0. 1582E-01 

0.2117E  00 

0.3902E  00 

0.7166E  00 

0.3699E-03 

0. 1687E-01 

-0.2940E-01 

0.3902E  00 

0.7649E  00 

0.1331E  01 

0. 1418E-02 

0.3100E-01 

-0.5389E-01 

0.7166E  00 

0.1331E  01 

0.2485E  01 

(C. 

2-1) 

Continuous -Time  Control  Rate  Weighting  Matrix 

[o-llllE  00  0.0 

R  —  I 

0.0  0 . 2778E-01 


(C.2-2) 


C-3 
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